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Abstract 

We report transport measurements on a quantum dot in a partly suspended carbon nan¬ 
otube. Electrostatic tuning allows us to modify and even switch “on” and “off” the coupling 
to the quantized stretching vibration across several charge states. The magnetic-field depen¬ 
dence indicates that only the two-electron spin-triplet excited state couples to the mechanical 
motion, indicating mechanical coupling to both the valley degree of freedom and the exchange 
interaction, in contrast to standard models. 
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Carbon nanotubes are found to be an ideal playground for nano-electromechanical systems (NEMS) 
since their high-quality, quantum-confined electronic states are accessible by transport spectro¬ 
scopic techniques and couple strongly to the excitations of different mechanical modes. The grow¬ 
ing interest in NEMS is fueled by the desire to accurately sense small masses and forces,^ ad¬ 
dress quantum-limited mechanical motion,!^ and integrate such functionality into complex hybrid 
devices,leading to new applications.® The central question is the strength of the coupling of elec¬ 
tronic states to the vibrational modes. Whereas molecular junctions display such modes also in 
electrically gated transport measurements,®® carbon-nanotube (CNT) quantum dots allow for a 
much more viable fabrication, higher mechanical Q-factors, and better tuneability as NEMS.®^ 
Also, the coupling to the bending mode can be combined^^®^ with the spin-orbit (SO) interac- 
tionChUllby making use of the recently demonstrated^ curvature-induced SO-coupling in CNTs.^ 
Whereas the frequency of the vibrational modes has been demonstrated to be tuneable,^^!^ an¬ 
other desirable feature is the ability to switch “on” and “off” the electron-vibration coupling in the 
same device, e.g., in envisioned quantum-information processing schemes.^^^^This is also helpful 
for fundamental studies of systems in which mechanical motion is combined with other degrees of 
freedom, e.g, the spin^^and the valley 

Recently, switchable coupling to a classical flexural mode of a CNT has been demonstrated.^ 
In this letter, we present a CNT quantum dot NEMS with a coupling of the electronic states to a 
longitudinal stretching vibration of about 200 GHz that can be turned “on” and “off”. We illustrate 
the advantage of this by transport measurements in the two-electron quantum-dot regime and find 
that the well-known Anderson-Holstein scenario breaks down in an unexpected way: Different 
spin states exhibit different coupling strengths to the vibrational mode. 

In Figure [T^ we show a schematic of a typical suspended CNT quantum-dot device whose 
scanning electron microscope image is shown in Figure [T]?. The CNT is electrically and mechan¬ 
ically connected to both source (s) and drain (d) contacts where the central electrode acts as a 
suspended, doubly clamped top gate (tg). The quantum dot is formed in the small band gap CNT 
by the electrostatic potentials of the top and back gate (bg), see Figures [T^ and |^, allowing for 
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Figure 1: Carbon-nanotube quantum dot characterization, (a) Schematic illustration of the cross- 
section of a partly suspended CNT connected to source and drain electrodes (5nm Cr, 50 nm 
Au). While the back-gate shifts the entire potential of the whole structure, the top-gate bridge 
overlaps with a 200 nm part of the CNT by an oxidized Cr layer, see Supporting Information, 
(b) Scanning electron-microscopy image of a partly suspended CNT sample. The CNT is just 
visible and indicated by red arrows, (c) Source-drain current through the quantum dot at zero 
magnetic field as function of the bias (Tsd) and top-gate voltage (Vtg), adjusting the back-gate 
voltage Vbg simultaneously to keep the average chemical potential in the leads constant: V^g = 
4.35 V — 0.7 X Vtg- (d) dl/dV^^ at zero magnetic field centered around the 1 2 single-electron 

tunneling regime for back-gate voltage Vbg = 4.25 V in the hole regime. We count the number of 
electrons relative to the last filled conduction band shell of the CNT as usual. The diagonal dashed 
lines marked A-F correspond to transitions between the N = 1 and 2 electron quantum-dot states, 
(e) Energy diagram of the one- and two-electron quantum dot states involving the first orbital shell 
and the corresponding orbital fillings discussed in the text. These states are responsible for the 
transitions A (Z)_ -H- 5-), B (D- Tfi) and C (D- ^o) which are the most relevant ones for 

the present discussion, see also Supporting Information. For panel (d) we extract A = 0.8 meV 
and 7=1.2 meV. (f) Measured magnetic field dependence of the electronic excitation lines along 
the horizontal line in panel (d) at Fsd = 6 mV. This way of plotting (see Supporting Information) 
focuses the attention on the important triplet states by making the transition B into the lowest triplet 
appear as a vertical line. 
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electrostatic control of the size of the quantum dot in the range of 250 - 350 nm, see Supporting 
Information. By changing the gate voltages we can modify the position and size of the dot with re¬ 
spect to the suspended vibrating region of the CNT, which is a crucial part of our experiment. The 
high quality of our CNT sample is revealed by the observation of well-resolved, multiple four-fold 
shell-structure of the electronic states in the stability diagram in Figure [T]: measured at zero mag¬ 
netic field and at a base temperature of 1.6 K. This shell-structure stems from the combined spin 
and valley degeneracies in clean CNTs,l2^^ and enables a first characterization of the electronic 
properties by the Coulomb and confinement energies. Importantly, the resulting estimates show 
that the quantum dot formed in the CNT is comparable to or even larger in size than the top-gate, 
see Supporting Information. 

The key advantage of our device, in contrast to previous ones, is that we can first obtain de¬ 
tailed information about the electronic spectrum by measuring the differential conductance in a 
gate voltage regime without signatures of vibrational excitations. For example, in the spectrum 
shown in Figure [T]3 the low-energy excitations indicated by dashed black lines can be assigned to 
transitions between states with electron number N = 1 and 2, respectively. These are indicated 
in the schematic in Figurewhich shows for N = 1 two spin doublets denoted D- and D+, ob¬ 
tained by filling the (anti)bonding orbitals |±) = (|^) ± |^'))/v^ of the K and K' valleys with 
one electron, which are split in energy by 2A due to the valley-mixing A. For N = 2 we have 
spin-singlets S- and 5+ (latter not shown) completely filling one of these orbitals, and a singlet So 
and a triplet Tq in which two different orbitals are filled. Here the labels of the many-body states 
S, D, T indicate the spin multiplicities (singlet, doublet, triplet), whereas the subscripts indicate 
the relevant orbital polarizations. In the transport data of Figure we identify a ground singlet 
(A-), an excited triplet (Tq) and another singlet (Aq), split by the exchange energy J. The measured 
magnetic field transport spectroscopy in Figure confirms this assignment: the slope of the lines 
A and C for transitions to A_ and Aq, respectively, differs by the Zeeman spin splitting from the 
slope of line B for the transition to the triplet Tq. We note that for these parameters the singlet A+ 
is the highest in energy in Figure [T^. It is not shown there nor discussed further below because this 
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Figure 2: Switching the coupling to the vibration “on” and“off”. (a) Schematic illustration of the 
quantum-dot tuning into a region with a longitudinal stretching mode (LSM). (b) Top- and back- 
gate voltage stability diagram recorded for 14(i = 1 mV. The lines marked (1) [Vbg = 4.41 V — 0.5 x 
Vtg] and (2) [Vt,g = 4.35 V — 1.0 x Vtg] indicate different regimes of electron-vibration coupling, (c) 
dZ/dVsd measured along line (1) in panel (b) showing no effects of vibrations, (d) Measurement 
of dl/dVsd along line (2) in panel (b), where significant electron-vibration coupling is observed: 
the arrows indicate the vibrational sidebands introduced. Electronic lines A and B from panel (c) 
can still be identified, but C is commensurate with a vibrational sideband of B. [Note that the 
same happens in the calculations in panel (g).] The blue markers indicate the end-points of the 
line (not shown for clarity) along which the measurements in Figure are taken, (e) Si mi lar 
measurement as in panel (d) but for a different relation of gate voltages (Vtg = 4.45 V — 0.95 x Vbg) 
showing vibrational excitations (arrows) with different gate-voltage slope, both in magnitude and 
sign. (f),(g) Calculated dl/dV^d corresponding to panel (c) and (d), respectively, see text. The 
overall conductance magnitude is adjusted through the coupling T, taking T = 0.7 K. (h) Linear 
fit of the vibrational excitations: the upper panel fits data from Figure]^ and the lower panel from 
Figure]^, both confirming a harmonic spectrum with AEvib = ^tt) = 0.8±0.1 meV corresponding 
to 193 ±24 GHz. 
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state does not influence the measured transport in the considered regime.l^Our calculations below 
do, however, include the state S+ and confirm that it has negligible influence. 

By independently tuning the top- and back gate voltages we can change the electrostatic con¬ 
finement of the quantum dot and thereby effectively operate a single quantum dot system which 
can be made sensitive to the vibrating part of the CNT, as illustrated in Figure |^. The result¬ 
ing electronic stability diagram in Figure |^, showing nearly parallel lines, indicates that we can 
independently fix the electron number in the dot while modifying its shape, dimensions and posi¬ 
tion. When measuring the Coulomb diamonds along the lines indicated in Figure]^ one expects, 
electronically speaking, no qualitative difference. Indeed, along the initial working line marked 
as (1) in Figure]^, the measurement in Figure]^ shows no indications of vibrations. However, 
when tuning to the working line (2), the excitation spectrum, shown in Figure |^, changes in a 
way that cannot be explained by a modification of the size-quantization energy on the quantum 
dot: for several subsequent charge states a dense spectrum of discrete excitation peaks appears, 
equally spaced by h(0 = 0.8 ±0.1 meV as Figure]^ shows. This is the case across the entire elec¬ 
tronic shell that we measure, see Supporting Information. The spacing lies in the range expected 
for the high frequency of the longitudinal stretching mode (LSM) of the suspended parts of the 
CNT (length ~ 65 nm as in previous studies^^^^^. Furthermore, the predominance of the excitation 
lines with negative slope indicates that the quantum dot couples to only one of the two suspended 
parts .1221201 demonstrate that by tuning to a different voltage regime we are able 

to make the other vibrating part dominate. Our system thus displays electrostatically tuneable 
electron-vibration coupling. 

To illustrate how the switchable coupling to a quantized vibration can be exploited, we now 
focus on measurements for the W = 1 -H- 2 electron regime in Figures ^ - 1^. The corresponding 
calculations shown in Figures |^ -are based on a quantum-dot model including the electronic 
states identified before in Figures [TJi - [T^ and coupling to a single vibrational mode. This model 
will be discussed in detail below, once we have presented all experimental data. Apart from this, 
the non-equilibrium transport is obtained from standard master equations (see Supporting Informa- 
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Figure 3: Magnetic field evolution of the vibrational sidebands, (a) dl/dVsd as function of top-gate 
voltage Hg along the line (not shown) connecting the blue markers Figure]^. Each dl/dV^d line 
has been measured at a different magnetic field B, which has been tuned from 5 = 0 to 9 T in steps 
of 300 mT. The dl/dV^d lines are offset vertically by 0.015 /h for clarity. Dashed lines indicate 
the singlet (sloped) and triplet (vertical) transitions, (b) Bias-voltage spacings (AFsd) of peaks A-E 
in panel (a) as function of the magnetic field including additional intermediate line traces that are 
not shown there. The dashed line for the peak-spacing A corresponds to g-factor 2. The data are 
offset vertically by 1.5 mV for clarity, (c) Calculated dl/dV^d evolution corresponding to panel (a). 
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tion) which incorporate single-electron tunneling into both orbitals of the shell (with asymmetry 
parameter k) from both electrodes (with junction asymmetry parameter 7 ). The electronic and 
vibrational states are assumed to relax with a phenomenological rate which exceeds the tunneling 
relaxation rates, taken for simplicity to be proportional to the energy change E in the transition: 
Trel (E) = rx {E/ 0.2 meV). The overall tunneling rate F merely sets the scale of the current and 
is irrelevant to the relative magnitude of the different excitations which is of interest here. 

To experimentally identify the electronic states to which the vibrational excitations belong, we 
have investigated how the differential conductance measured along the line (not shown) connecting 
the blue markers in Figure evolves with a magnetic field B applied perpendicular to the CNT. 
The dominating features in Figure are the vibrational sidebands of the lowest of the triplet 
excitations Tq which in this presentation of the data appear as vertical lines. Strikingly, the ground 
state singlet 5_ evolving with a slope has no vibrational sidebands as demonstrated by fits of 
the difference of the peak position in Figure |^. This can not be explained by an Anderson- 
Holstein model where all electronic states with the same charge couple equally to the vibration, 
see Supporting Information for explicit attempts. 

Instead, in our modeling we must account for state-dependent Franck-Condon shifts resulting 
in the vibrational potentials plotted in Figure]^ To arrive at this, we start from a model accounting 
for the observed set of accessible!^ many-body transport states, which is restricted by Coulomb 
blockade and bias voltage of a few mV to those shown in Figure [T^ with electron numbers N = 1 
and N = 2 and a single electronic K-K' shell: 

//ei = eA^ + A£ Tj^4^d,cT-/S+-S_. (1) 

T=± (7 

Here e is mean level position controlled by Vtg, A is the subband or valley-mixing term and J 
is the exchange coupling between the spins in the two orbitals T = ± with spin-operators = 
5 Lcr,cr'<^T <7 0'cT,<7'^TCT' spin indices, a is the usual vector of Pauli-matrices, and 4<7 cre¬ 

ates a spin-(T electron in orbital t]. To obtain a result as plotted in Figures]^ -|^ we first introduced 
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a Holstein coupling by allowing the level position e to depend on Q, the dimensionless vibration 
coordinate normalized to the zero-point motion: we thus formally replace £—>■£ + 

This results in the commonly assumed uniform vibration coupling with strength A-e to all elec¬ 
tronic states with the same charge N, which is not observed here. The required state dependent 
electron-vibration coupling is obtained by additionally accounting for a dependence of the other 
parameters on the vibration coordinate, i.e., we formally replace A ^ A-|- \/2h(oX/^Q, where Aa 
is a dimensionless one-electron valley-vibration coupling, and J ^ J + '/ITkdXjQ, where Ay is a 
dimensionless coupling of the vibration to the two-electron exchange. Here many-body physics 
comes in: when going from the singlet 5_ ground state to the triplet Tq, the Pauli principle forces 
the two electrons into different orbitals which can couple differently to the vibrational mode (dif¬ 
ference quantified by Aa). However, the coupling Ay is important as well: when allowing only for 
Aa, the effective electronic excitation spectrum for fixed charge N (relative to which the vibration 
excitations are “counted”) becomes dependent on the vibrational couplings (polaronic renormal¬ 
ization). That experimentally no significant shift of the electronic excitations is found when turning 
“on” the couplings to the vibration requires the couplings Ay and Aa to be comparable in magni¬ 
tude but opposite in sign. This results in an enhanced coupling of the triplet Tq over S- while 
the polaronic shifts that they induce cancel out, keeping the effective electronic excitations fixed. 
This thus leaves one free parameter, their magnitude, which controls the degree of state-specific 
coupling, which we adjust to the experiment. Together this suffices to obtain results such as Fig¬ 
ures B-a that reproduce the main zero-field observations of Figures]^ -|^. When the vibration 
couplings are “off’ in Figurewe estimate from Figure]^ the parameter values A = 0.8 meV, 
/ = 1.5 meV, (similar to those in Figure [T]l) and use 7 = 0 . 0 , K = —0.3. When the vibration 
couplings are “on” in Figure we use the same values for J and A but nonzero vibration cou¬ 
plings Ag = 0.28, Aa = 0.32, Ay = — 0.22 and frequency hco = 0.85 meV and we adjusted the 
asymmetries 7 = —0.5, K = 0.3. Despite the fact that there are several parameters, the experiment 
imposes strong restrictions, in particular, regarding the choice of vibrational couplings, excluding 
a simple Holstein mechanism (Aa = Ay = 0), see Supporting Information. We arrive at the three 
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electron-vibration couplings by imposing three experimental constraints after expressing the effec¬ 
tive couplings of the electronic states in terms of Ae, Aa, and Ay: (i) the observed Tq-S- splitting 
and (ii) So - Tq splitting (commensurate with 2h(0) should match energy expressions that depend 
on the vibrational couplings (polaron shift) and (iii) the vibrational-coupling of Tq is adjusted to 
numerically reproduce the observed number of triplet vibrational sidebands. We note that in Fig¬ 
ures]^ -|^, the higher vibrational sidebands become more intense at high bias. As expected, this 
is not captured by our model since this may involve excitations beyond the lowest two electronic 
orbitals and energy-dependence of the tunnel barrier, neither of which we include here. We have 
focused instead on the nontrivial interplay of vibrational and spin-excitations for A = 1 and N = 2 
in the lowest sidebands. 

The resulting physical mechanism is illustrated in Figure when starting out from state D- 
and adding a second electron to the lowest orbital the lowest singlet state S- experiences only 
a small horizontal shift of the vibrational potential minimum (both electrons in orbital |—) have 
their coupling weakened by Aa and there is no spin and therefore no exchange modification of the 
coupling by Ay). However, when adding the electron to the excited orbital, the coupling is not 
only enhanced by Aa, but also by a negative Ay when a spinfull triplet Tq is formed. This results 
in a large Franck-Condon shift of the potential minimum of Tq in Figure The above horizontal 
shifts of the potential minima translate into suppressed vibrational sidebands for the singlet S- 
and a pronounced series of sidebands for the triplet Tq, respectively (Franck-Condon effect). The 
presence of the further electronic states and their quantized vibrational states in Figure all of 
which are included in our transport calculations, do not alter the above simple picture: Whereas 
the excited singlet S+ does not couple to transport,!^ the role of Sq cannot be ascertained at zero 
magnetic field because it is commensurate (within the line broadening) with one of the vibrational 
sidebands of Tq. 

The field evolution in Figure]^, calculated by adding a Zeeman term to equation (??), repro¬ 
duces the main observation of Figurej^, namely, that the triplet maintains its vibrational sidebands 
(vertical) but the ground singlet S- (sloped) does not. However, to obtain this agreement with the 
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Figure 4: Vibrational potential energy surfaces underlying the quantum states of the vibrating 
CNT quantum dot included in the transport calculations using the same parameters as in Figures 
- 1^ and [^. The excited triplet Tq and ground singlet 5_ have significantly different couplings 
to the vibration, i.e., shifts of their potential minima relative to that of the one-electron ground 
state D- strongly differ. Due to the weak SO coupling several avoided crossings can be seen. The 
most important anticrossing is that of the Tq (blue) and 5_ (red) potential energy surfaces. This 
can be understood directly from the SO operator as written in the text: it "flips" both the orbital 
(t) and spin index (cr) of an electron. In Figure [T^ this implies that for V = 2 the red spin up 
in the higher orbital is flipped into the blue spin down in the lower orbital (this represents a flip 
from Tq to 5_). The resulting admixture of To-components (blue) to S- (red) causes the latter to 
remain visible in the transport in Figure]^ with increasing the magnetic field when the tunneling 
becomes spin-selective due to the CNT leads. The remaining SO anticrossings are discussed in the 
Supporting Information, which for our parameters have negligible impact on transport. 
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measurements we are forced to further extend the above model. First, both the excited singlet 
(^o) as well as the Zeeman split-off states of the triplet (Tq) do not appear in the measurements. 
This we attribute to the fact that the source and drain leads of the quantum dot are not formed 
by metallic contacts but by small pieces of suspended CNT. Zeeman splitting of discrete states in 
these CNT contacts may lead to spin-filtering which turns on with the magnetic field, developing 
full strength at a few Tesla where gfi^B ~ We phenomenologically account for this by a spin- 
dependence in tunneling to / from the electrodes which depends on B: ^{B) = tanh(gjUB5/2^BT). 
Second, when only including this spin-filtering in the model, it suppresses the singlet groundstate 
S- (without vibrational bands) which we do experimentally observe as excitation A in Figures]^ 
-a. However, when even a small spin-orbit (SO) coupling is included, the singlet 5_ reappears 
(borrowing intensity from the triplet Tq, cf. also Figure]^, but, importantly, without reinstating the 
unobserved 5o and the Zeeman split-off states of Tq and their vibrational sidebands. This produces 
the observed intensity pattern, which is impossible to achieve with simple commonly used models, 
see Supporting Information. Here, the spin-orbit coupling is included by adding to equation (??) 
a term Hso = ^soY.a,T<^lad-T-a with Aso = 0.1 meV which allows both the spin G and orbital 
index T to be flipped in the schematic Figure [T]l, thereby coupling in particular Tq to 5_, lending it 
intensity. Figures]^ -|^ and|^ are based on the inclusion of all these effects. However, we empha¬ 
size, that in the latter figure spin-filtering and spin-orbit coupling are needed exclusively to explain 
the missing Zeeman lines, but do not lead to a qualitative change of the state-dependent coupling at 
= 0 in Figure]^, which is our main finding. The Supporting Information explores the influence 
of the various parameters, confirming the necessity of including them. The key advantage of our 
tuneable setup is that we are able to first identify excitation A and B as relating to electronic sin¬ 
glet S- and triplet Tq, respectively, and subsequently allowing us to study the vibrational sidebands 
C-E. 

In conclusion, we have demonstrated switchable coupling of a quantized vibration of a carbon 
nanotube to its quantized electronic states. Using this advance we explored the two-electron regime 
- including the magnetic field dependence - and found indications of state-dependent vibrational 
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transport sidebands not described by standard models. We showed that the interplay of intrinsic 
effects on the carbon nanotube (Coulomb blockade, valley-index, spin-exchange) and experimen¬ 
tal details (junction, orbital, and spin asymmetries) can explain the observations. This, however, 
includes vibrational couplings that involve internal spin- and valley-degrees of freedom, bringing 
spin- and valley-tronics physics within range of NEMS. 
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In this Supplementary Information we provide a detailed description of the experimental meth¬ 
ods and additional measurements (Sec. I) and a precise formulation of the model and calculations 
reported in the main text (Sec. II). In both sections, we provide an extensive discussion of claims 
and results of the main article. Within the Supporting Information references are numbered as, 
e.g., equation (S-1) and Figure S-1, whereas regular numbers, e.g., equation (1) and Figure 1, refer 
to the main article. 
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(a) 


FIG. S-1. Fabrication work flow described in the text. 

I. EXPERIMENTAL METHODS 
A. Fabrication 

Devices were fabricated in a similar fashion as in Ref. [1] as outlined in Fig. S-1. The starting 
point is a highly doped silicon wafer covered by 290 nm silicon oxide (step a). Next, following the 
recipe described in Ref. [2] ferritin catalyst nanoparticles are dispersed on the substrate (step b) 
from which carbon nanotubes are grown by means of chemical vapor deposition (step c). For the 
subsequent selection and localization of carbon nanotubes marker structures are evaporated in an 
electron beam (e-beam) lithography step (step d,e). In a second e-beam lithography step metallic 
electrodes and gate structures are deposited in a single evaporation (5nm Cr, 50nm An) on selected 
carbon nanotubes (step f-h). Finally, diluted hydrofluoric acid (1% for 6min) is used to etch the 
silicon oxide followed by critical point drying (step i). It is crucial that the central electrode is 
completely underetched so that the chromium layer oxidizes when exposed to environmental air 
to form the top-gate oxide. 




B. Configuration of the quantum dot 

In Fig. S-2 we show the back-gate voltage (Rbg) characteristics of the investigated carbon nan¬ 
otube (CNT) at two different temperatures. The CNT is slightly p-doped and the charge neutrality 
point is found to be around 7 V. At gate voltages close to the charge neutrality point the current 
Jsd is strongly suppressed, even for elevated temperatures T = 50 K. This is the typical character¬ 
istics of a CNT with a small semiconducting gap separating p- and n-conducting regions. For low 
temperatures the current within the semiconducting band gap, apart from a few small resonances, 
is pinched off. The band gap extends over a back-gate voltage range of Al^g = 1.5 V. Multiplied 
with the absolute lever arm of the back gate o^bg = 0.035 the energy splitting between valence 
and conduction band is AEgap ~ 50meV. At T = 1.6 K we observe a number of reproducible 
resonances. These resonances are attributed to Coulomb blockade effects. The irregularity of 
the resonances indicates the existence of multiple quantum dots and a non-monotonic CNT band 
structure along the nanotube axis. 

In Fig. S-3a we illustrate the electrostatic formation of the main quantum dot. The Fermi level 
along the entire CNT can be tuned with respect to valence and conduction band by an applied 
back-gate voltage Vbg- Owing to the geometry of the device (see Figures la - lb of the main text) 
different electrostatic potentials act on different parts of the nanotube depending on the proximity 
to the local top gate: The electrons just below the top gate experience electrostatic screening 
of the back-gate voltage because the CNT is separated only by a few nanometers of oxide from 
the top-gate electrode while a combination of top gate and back-gate voltage is acting on the 
suspended parts of the CNT. Therefore, we expect a bending of the CNT band structure along its 
axis. 

A device schematics together with an indication of the location and size of the quantum dot 
with respect to the metallic leads and the top gate is given in Fig. S-3b. As deduced below, the 
quantum dot is formed close to the top gate and its lateral size is on the order of the top-gate 
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FIG. S-2. Source-drain current Isd through the CNT as a function of back gate voltage Vbg at small bias 
voltage Vsd = 0.3 mV at two different temperatures. 

(a) (b) 




Back Gate 




FIG. S-3. Formation scenario of the quantum dot. (a) Electrostatic bending of the GNT bands in close 
vicinity to the top gate induced by top and back gate voltages, (b) The upper illustration shows a 
schematic of the device. The semi-suspended GNT is connected to Gr/Au leads and a central top gate. In 
the suspended region longitudinal stretching modes (LSM) of vibration are indicated by the blue spirals 
surrounding the GNT, while the red region below the top gate symbolizes the quantum dot. The lower 
schematics describes the tunneling through the quantum dot. Electrostatically induced tunnel barriers 
separate the GNT quantum dot electrically from the GNT leads. 


width of 200 nm. These parameters suggest, that the leads of the QD are not the metallic contacts 
but the CNT itself and therefore the size of the QD can be tuned exclusively by the applied gate 
voltages. 

In order to determine the electronic size of the quantum dot we analyse the addition energies 
needed to add the first and the second electron on the investigated electronic shell of the QD, 
respectively. The addition energy is defined as the change in electrochemical potential A/i 7 v when 
adding the (A + 1) charge to a quantum dot containing already N charges^’^. It can be related 
to the charging energy the quantum energy-level separation (3^ the valley degeneracy splitting 
2A and the exchange interaction J via 

A/io = ^c+;5-2A (S-1) 

A/ii = Ec- (S-2) 

Additionally, we can infer J and A from the energies of the first two excited states relative to 
the two-electron ground state, Etq — Es_ = 2A — J/4 and Esq — Es_ = 2A -f 3J/4, respectively 
[cf. equation (S-14), (S-45), and (S-46) below)]. From Fig. 2c of the main article [reproduced in 
Fig. S-6b below] we extract A/tq = 9.0 ± 0.3 meV, A/ii = 3.8 ± 0.4 meV, /ii = 1.2 ± 0.3 meV and 
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FIG. S-4. (a) Charge stability map in the Vtg-Vbg-plane. The differential conductance is plotted for 

zero d.c. bias voltage Vsd applied between source and drain. Resonances with a slope marked by the 
blue dashed lines are attributed to resonances off the central quantum dot. The Coulomb diamonds in 
Fig. Ic of the main article are measured along the white continuous line and the Coulomb peaks in (b) 
are measured along the white dashed line (Vbg = 4.26 V — 0.6 x Vtg)- 


/i 2 = 2.7 ± 0.3 meV. This yields 

[3 = 6.8 meV, J = 1.5 meV and A = 0.8 meV. (S-3) 

Due to large level spacing (3 relative to J and A, the analysis of the transport spectrum involving 
the first few vibrational excitations - at the focus of the main article - only requires considering 
electron fillings of the first orbital shell (i.e., of the orbitals labeled r = ± in the main article). A 
more detailed fitting of the energy positions consistent with the above and including the vibrations 
is given in Sec. IIC 1. The expression for the quantization-induced level spacing /3 = with 

v-p = S.l X 10^ m/s allows to determine the quantum dot length L = 245 ± 25 nm, corresponding 
approximately to the length below the top gate. 

For the determination of the position of the quantum dot we have evaluated the respective lever 
arms of both top gate and back-gate electrodes as well as the relative lever arms of the source and 
drain leads. The analysis of the slopes of the edges of the Coulomb diamonds in Fig. Id of the 
main article gives the following lever arms: atg = 0.55 ± 0.04 and o^s — o^d = 0.01. 

Knowing the lever arm to the top gate, the lever arm to the back gate can be obtained from 
the charge-stability map when varying both top gate and back-gate voltages, as shown in Fig. S-4. 
This figure consists of lines with different slopes, which is a signature of a multiple quantum dot 
structure. The interesting features are the ones with the largest slope, marked by the blue dashed 
lines. The slope of these lines allows to extract the relative lever arm between top gate and back 
gate Qfrei = Q^tg/<^bg = 15.5 ± 1.9, resulting in o^bg = 0.035 ± 0.005. 
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FIG. S-5. Charge stability diagram in both gate voltages for finite d.c. bias voltage of Vsd = 1 mV. The 
differential conductance plot shows splitting of the lines belonging to the central quantum dot. 


The significantly larger lever arm to the top gate and the almost negligible difference between 
the lever arms of source and drain confirm our assumption that the quantum dot is located in 
the middle of the nanotube in close vicinity of the top gate and is as a result strongly screened 
from the back gate. Additional lines with smaller slope in the charge-stability diagram indicate 
the existence of additional quantum dots further away from the top gate. 

In Fig. S-4a we show a charge stability map as a function of both the top gate and the back-gate 
voltage. Strikingly, the Coulomb-peak excitation lines with the steepest slope, which correspond 
to the central QD, appear in groups of four, reflecting the twofold orbital-degenerate bandstructure 
of high-quality CNTs. Further confirmation of the fourfold periodicity is provided in Fig. S-4b 
where we plot a Coulomb-peak measurement as function of Ffg. These characteristics justify the 
treatment of the CNT QD as an effective few-electron system. 

The CNT-QD is thus connected to the metallic electrodes by short CNT leads. To prove that 
the CNT-leads do not electrically influence our measurements on the central QD significantly, we 
have measured the same diagram as in Fig. S-4a, but with a finite bias voltage Vsd = 1 mV. 
Fig. S-5 shows that only the lines corresponding to the central QD split up into two lines, while 
the lines corresponding to the lateral parts are not affected. This broadening into ground and 
excited states proves that the voltage drop occurs only at the tunnel barriers (marked in Fig. S- 
3b) between the CNT-QD and the CNT leads. The tunneling rates are estimated to be on the 
order of F = 27r x 10 GHz, which corresponds to roughly 300-400 mK, well below the experimental 
temperature T = 1.6 K. As a first approximation it thus makes sense to apply a standard master- 
equation description of single-electron tunneling transport, see Sec. II. 


C. Electrostatic control of the coupling to vibrational modes 

In this section we provide further experimental data on the tuning of the coupling of the QD to 
vibrational modes. In particular, in Fig. S-6 we provide an additional QD excitation spectrum to 
further illustrate our ability to continuously tune the electron-vibration coupling. Figures S-6a, 
S-6b and S-6d correspond to Figures 2b, 2c and 2d of the main article and show a pure electronic 
excitation spectrum and a vibrational excitation spectrum, respectively. In addition, we show in 
Fig. S-6c data for an intermediate state of the QD. 

In order to understand the switchable coupling of the QD to vibrational modes we compare the 
quantum dot size in the two extreme regimes. Unfortunately, when the vibrations are switched on 
we can not assign the energy of the second excited state So in the two electron regime (because 
it is degenerate with a vibrational sideband within the experimental line width). Nevertheless, 
we can give bounds for the dot size. With A/io = 9.3 meV, A/ii = 5.8 meV and /ii = 0.9 meV 
(see Fig. S-6d and also Fig. S-9b) we get P = 4.4 meV + J/4, which gives a larger quantum dot 
size in the vibrational regime for 0 < J < 9.6 meV. For realistic exchange energies this inequality 
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FIG. S-6. (a) Differential conductance in the Vtg-Vbg-plane for small d.c. bias voltage, (b)-(d) Coulomb 
diamond charge stability diagrams in the I4d-htg-plane for the effective charge states 0 — 4 in the quantum 
dot measured along the respective lines in panel (a). 


is fullfilled. If we choose, for example, J = 1.5 meV as in the regime without vibrations the 
quantum dot length exceeds 300 nm. The larger quantum dot size, when vibrations are switched 
on, indicates that the overlap between QD and suspended regions of the CNT is increased. This, 
in turn, results in coupling to vibrational modes. 


D. Temperature dependence of vibrational sidebands 


We also investigated temperature dependence of ground and excited states in the temperature 
regime where Coulomb blockade peaks still could be resolved along the lines of Ref. [1]. In Fig. S- 
7a we show Coulomb diamonds measured in the very same region of gate voltage (same electronic 
state) for different temperatures, 1.6 K, 2.5 K, 3.5 K and 5 K. As the temperature increases, 
the conductance peaks related to ground and excited states wash out, i.e., they broaden and 
the conductance maximum decreases. In Fig. S-7b we show the temperature dependence of the 
maximum conductance Gmax for the electronic triplet excited state Tq (blue circles) and its first 
vibrational replica, i.e. the emission side-band (red triangles). In the four panels in Fig. S-7a 
these excitations are marked by blue and red arrows, respectively. For the tunneling through the 
electronic triplet Tq excited state and its vibrational emission peak, we observe a Gmax ^ l/ksT 
dependence (blue and red curves) that one expects for the derivative of the Fermi distribution in 
the quantum Coulomb blockade regime^. 

Unlike in Figs. 2c-2d in Ref. [1] we do not observe any extra Yihiaiionol- absorption conductance 
peaks appearing with increasing temperature inside the Coulomb blockaded region, most likely 
because of weaker electron-vibron coupling strength / smaller tunneling rates. This prohibits 
further investigation of the temperature dependence. 
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FIG. S-7. (a) Coulomb diamond charge stability diagrams in the I4d-Ftg-plane around the N = 1-2 

transition for T = 1.6 K, as shown in Fig. 2d of the main article, and for three higher temperatures T — 
2.5 K, 3.5 K and 5 K. (b) Conductance peak value as function of temperature for electronic triplet state 
(blue circles) and the first vibrational replica (red triangles). The blue and red curves are guides to the 
eye oc l/Zc^T . 


E. Measurements in a magnetic field 

The motivation for plotting the data as done in Fig. 3a in the main article is explained in 
Sec. IIC 2. 


II. THEORETICAL MODELING 

In this section we describe in detail the employed model, the method used for calculating 
the differential conductance of the CNT quantum dot, and the resulting understanding of the 
transport measurements. The full Hamiltonian of the system under consideration reads as H = 
i7qd + i^tun + 77res5 whcrc i7qd dcscribcs the quantum dot states, including both their electronic 
and vibrational degrees of freedom, and i^tun is the tunnel coupling Hamiltonian between the dot 
and the reservoirs described by i^res- 


A. Model and eigenstates 

1. Electronic model 

a. Carbon-nanotuhe quantum dot We first set up an electronic model that accounts for the 
many-electron states observed in the experiment when the coupling to the CNT vibrations is 
switched “off”. In Sec. HA2 we then include the vibrations and their coupling to obtain our full 
CNT quantum-dot Hamiltonian i7qd. As mentioned at the end of Sec. IB we can restrict our 
attention to a single orbital shell. We account for a significant valley-mixing A > 0: 


AE(4a^^v+H.c.), 

a 


(S-4) 


where {dx'a) is a creation (annihilation) operator for an electron in valley K {K') with spin 
cr = {t, i} = {1, — 1} whose quantization axis is chosen along the direction of the applied magnetic 
field. Due to the large splitting A, it is reasonable to use a basis of bonding (r = —1) and 






























antibonding (r = +1) combinations (BA) of the K and K' valleys, i.e. dr} = 

In this basis, the electronic model Hamiltonian can be written as 

=eN- qi^bBS, + A ^ - JS+ • S_ + ^N{N - 1). (S-5) 

r 

In the first term, 5 = —atgVtg is the quantum-dot energy level, electrostatically controlled by the 
top-gate, and N = J^rcr is the occupation number operator in the dot. Note that in the 

experiment back- and top-gate are tuned simultaneously in a linearly dependent way and that 
Vtg is taken as the independent parameter. Here A = 1 — 4 counts the electrons that fill up the 
orbital shell that we consider and N = 0 corresponds to the ‘empty dot’ state |0) in which all lower 
shells are filled. The further number operators Nr = d\^dra and No- = count the 

number of electrons in orbital r and with spin a, respectively. The second term in equation (S-5) 
describes the spin Zeeman splitting due to magnetic field B applied perpendicular to the CNT axis 
and Sz = X]cr(^/^)^cr is the operator of the spin component along the field B. The third term in 
the electronic Hamiltonian (S-5) is a spin-exchange term with energy J. The spin operator S^- for 
electrons in orbital r = ± in the BA-basis is Sr = Z^crcr' ^cra'd^^dra 'where a is the vector of 
Pauli matrices. Finally, the last term is the charging energy Ec accounting for Coulomb repulsion. 
The inter- and intra-valley electronic repulsion energies are assumed to be the same as is typically 
observed in CNTs quantum-dot samples^. Since we will restrict ourselves to the analysis of the 
1 2 charge regions, both the charging energy Ec and the level energy p [cf. equation (S-1)] can 

be absorbed by a redefinition of the origin of top-gate voltage and therefore these do not need to 
be included henceforth. Thus setting Ec = 0 we obtain the Hamiltonian (1) of the main article, 
where we note that for brevity the Zeeman term was not written but only mentioned in the main 
article. 

A convenient many-particle basis for A = 1 and N = 2 charge sectors is constructed by creating 
spin-multiplets using these orbitals. In the following we will denote by \x^y) such electron fillings, 
where the first (second) slot x (y) corresponds to the ground orbital — (excited orbital +). By 
filling the empty dot state |0) with one electron, we obtain two spin-doublets denoted by D± in 
the main article, with states 


\DZ) = dl,\0) = \a,.), (S-6) 

\Dl)=dl,\0) = \.,a). (S-7) 

where • denotes an empty ± orbital, respectively. The next step is to add a further electron to 
these states. We therefore obtain two ‘localized’ singlet-fillings of the same orbital 

\S-)=dl/_^\0} = \if), (S-8) 

|5+)=4X^|0) = |.4t). (S-9) 

By filling the empty dot with two electrons in different orbitals and diagonalizing the electronic 
Hamiltonian of equation (S-5), we obtain a ‘delocalized’ singlet 


=-Vi-= —Vi— 

and a triplet Tq of states \T^) with spin-projections m = 0, ±1: 


1 ^0°) = 


dNU + 

V2 


't 


| 0 ) = 


ta) +1 'I'jt) 

V2 


|To+) = 4yqio) =lt,t), 

|To-) = <yq|0) =U4). 


(S-10) 


(S-11) 

(S-12) 

(S-13) 


The labels of the many-body states S, D, T indicate the spin multiplicities (singlet, doublet, 
triplet), whereas the subscripts indicate the relevant orbital polarizations (signature of difference 
of number of r = ± electrons, respectively), which is important here. 
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The above defined states are exact many-particle eigenstates of the model (S-5) whose corre¬ 
sponding eigenenergies are 

r Es, = 3J/4 

N = l: = tA — /2 N = 2: < Et^ = — (S-14) 

[ Es^ = 2rA 


For the moment we consider zero magnetic field 5 = 0 and can simplify the discussion by omitting 
the spin projection indices a in the doublets and m in the triplet. In Sec. IIC we will return to 
this notation when discussing the effect of a nonzero magnetic field. Clearly, in the N = 1 charge 
sector the ground state is the bonding doublet D_ while for A = 2 the ground state is given 
by the localized singlet S- in the expected regime of weak exchange energy relative to the valley 
mixing, J < 8A. 

b. Transport model The source (s) and drain (d) leads are described as macroscopic reservoirs 
of noninteracting electrons through the Hamiltonian 

Hres ~ ^ T Tr)^rk<j^rk(j i (S-15) 

rka 


where {crko) creates (annihilates) an electron in lead r = {s, d} with spin a = {t, 1} and state 
index k. The eigenenergies of the leads are uniformly shifted by the bias voltage Vsd such that 
the electro-chemical potentials read (ir = ±I4d/2 for r = {s,d}, respectively. These reservoirs are 
assumed to be independently at equilibrium, characterized by a temperature T. We note that the 
matter of interactions in the CNT leads is a subtle one, but in the experiment we see no particular 
effect indicating their importance. Rather, the fact that we have quantized (yet broadened) states 
in the CNT leads seems to be important, resulting in an effective spin-dependence of the tunneling 
rates, see below. 

The coupling between the dot and the leads is determined by the tunnel Hamiltonian 

fftun = ^ KAaCrka + H.C. , (S-16) 

rkra 

with the tunnel amplitudes assumed to be junction- (r), orbital- (r) and spin-dependent (a). 
Since the tunnel rates required below have the form 27rp(t^^)^ (by Fermi’s Golden Rule) these 
dependencies are modeled using three asymmetry parameters k, 7 and C, respectively: 

C = + + (S-I7) 

where r = ± corresponds to r = s/d. Here E characterizes the tunneling through each junction r 
through the overall rates F’^ = 27Tp{E)‘^ and p is the density of states in the respective electrode. 
We let F = F^ + F^ eharaeterize the overall seale of the rates whieh merely sets the magnitude 
of the eurrent and is irrelevant to the relative strengths of the different excitations which are of 
interest here. The latter are controlled by the quantum dot electron-vibrational states and the 
parameters 7 , /^, and 

• In equation (S-17) we include the usual junetion asymmetry 7 = (t^ — t^)/{E -h t^) for the 
tunnel eoupling to the souree and drain leads through the last factor. Keeping F = F^ -h F"^ 
fixed, a little algebra shows that E is given by the second equation in equation (S-17). This 
asymmetry is relevant for modeling the measured differential conductance, which manifests 
some asymmetry between the intensities of positively and negatively sloped lines. We note 
that any r-dependence in the density of states in the leads that we ignored above can be 
absorbed into 7 . 

• The first factor in equation (S-17) captures an orbital asymmetry k = “^-cr)/(^+cr+^-cr) 

which is the same for all a and r. This derives from the linear combinations of the valley- 
dependent tunnel amplitudes, i.e. = {t^^ + E_^)/^/2 and = {t^^ — E_^^l\f2. We 
notice that for strictly symmetric couplings of the two orbitals, k — E^'al'^Ka ~ tunnel 

Hamiltonian is symmetric with respect to the interchange of the K and K' valleys and cannot 
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induce transitions from the anti-symmetric N = 1 ground state D_ into the symmetric state 
N = 2 ground state S-. This would cause a strong suppression of low-bias transport up to 
a voltage where the lowest excitation for either N = 1 or N = 2 becomes accessible. This is 
not observed in the experiment and indicates that a definite orbital asymmetry is present, 
i.e., n ^ ^1. (A similar problem arises for n = —1.) 

• Finally, we introduced in addition a possible spin-dependence in the tunnel amplitudes 
through the parameter ( = {t^^ — t^^)/{t^^ -h t^^) which is the same for all r and r. Below 
we turn on ( only in a nonzero applied magnetic field B [cf. equation (S-59)]. As mentioned 
in the main article, this models a relevant aspect of the experiment, related to the fact that 
we have CNT leads, which we will discuss in Sec. IIC when calculating the magnetic field 
evolution of the vibrational sideband lines observed in Fig. 3a of the main article. 

From the above Hamiltonian we calculate the tunnel matrix elements (TMEs) 

T:il = Y,tlMdUb), (S-18) 

r 

with the shorthand d^^ = d\^^dra for = ±1. In Sec. IIB we derive the explicit TMEs for the 
considered model after having considered the effect of the vibration, to which we turn now. 


2. Coupling to the vibration - beyond the Anders on-Holstein model 


The electronic model accounting for the many-electron states observed in the experiment is now 
extended to deal with the case where the coupling to the vibrational stretching mode of the nan¬ 
otube is switched “on”. As discussed in the main article, the differential conductance of Eig. 2d 
together with its magnetic field evolution in Eig. 3a shows several vibrational sidebands associated 
to the triplet state, but none for the ground singlet. This strong state-dependence in the coupling 
to the vibrational mode forces us to consider three different types of electron-vibration couplings 
which arise from the assumption of a linear dependence on the (dimensionless) mechanical dis¬ 
placement Q of the nanotube in the electronic parameters R = { 5 , A, J} of equation (S-5). In all 
cases, we assume R{Q) = R-\- \RhLJ\/2Q. We therefore have, in addition to the standard Holstein 
coupling Ae to the number of particles A, a coupling Aa which depends on the valley-mixing and 
a vibration-exchange coupling Aj. By plugging these into equation (S-5) we arrive to the full 
Hamiltonian of the CNT quantum dot 


hjjj r 


Hqd — He\ T ^ 


p2 + (Q + ^/2Kf 


hujA^ 


(S-19) 


which can be seen as a shifted quantum harmonic oscillator. The above electron-vibration cou¬ 
plings enter through the following operator 


A = X,N + XAYi - AjS+ • S_ , 

T(J 


(S-20) 


which shifts the harmonic potentials associated to each electronic state (horizontal shift) and it 
also introduces a polaronic shift in the energy (vertical shift). 

a. Adiabatic potentials To aid the intuition we consider the adiabatic potentials for this prob¬ 
lem, obtained by treating Q as a classical variable. These potentials were plotted in the main article 
in Eig. 4. The adiabatic potentials associated to the electronic states |e), where e = D±^ S±^ Sq^ Tq, 
can be characterized by polaron shift of its potential minimum Ae = (e|A|e), i.e. 

Ad^ = Ag + rA a , As^ = 2\^ + 2rA a , (S-2I) 

= 2A£ + -Aj , Ato = 2A£ - -Aj . (S-22) 

Roughly speaking, the magnitude of these shifts determine whether none, several or many vi¬ 
brational sidebands will appear (see below). These expressions reveal that the polaronic shifts 
corresponding to the various electronic states can indeed be different once one abandons the 
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simplifying assumption of the Anderson-Holstein model, that the electronic excitations are not 
affected by a distortion (the vibration coordinate Q)\ 

(i) Both the valley-vibration coupling A a and the exchange-vibration coupling A j distinguish mul- 
tiplets with zero {Sq^Tq) and maximal orbital polarization {Dr,Sr)- 

(ii) Also, the exchange-vibration coupling allows for a fine tuning of the coupling to the vibration 
for the So and Tq states. As we will show in Sec. IIC, the horizontal shifts of the ground singlet 
S- and the triplet Tq can be tuned so they are qualitatively different, giving rise to contrasting 
Franck-Condon factors. 

Although the intuition is useful, in our calculations we account for the exact many-body eigen¬ 
states of the quantum-dot Hamiltonian in equation (S-19). These are obtained as a tensor product 
of the electronic states |e) and the quantum vibrational states \n). These electron-vibration eigen¬ 
states |e, z/) thus yield the following eigenenergies for e = D±, S±, So, Tq: 

Ee,v = (e, i^|i?qd|e, v)= Ee - fujAl + ^ • {S-23) 

The horizontal shift in equation (S-19), although not present in the above energies, enters as the 
equilibrium position for the vibrational states, and thereby it plays a crucial role in the transition 
amplitudes when calculating the matrix elements of the tunnel Hamiltonian, as we will show in 
Sec. HB. 

Finally, we note that in the limit Aa = Aj = 0 the model reduces the standard Anderson- 
Holstein model. Here, the polaronic energy shifts can be absorbed into the gate voltage since the 
electron-vibration coupling for transitions involving states differing by one electron is, in all cases, 
the same. 

b. State-dependent Franck-Condon shifts - the role of Xa cind Xj The tunnel amplitudes and 
thus the intensities of the lines in the differential conductance depend on the Franck-Condon 
amplitudes that can be calculated from the relative shifts between the harmonic potentials, i.e. 
Ax^y = Ax — Ay where x, y are labels of the many-particle states differing by one electron. Of 
particular relevance to the experiment are the two relative oscillator shifts 

As_,d- = Ae - Aa , (S-24) 

^To,D_ = Ae + Aa - -Aj , (S-25) 

for transitions from the ground doublet D- to the ground singlet S- and to the triplet Tq, 
respectively. Since Aa is present in both transitions but with different sign, it can induce different 
coupling strengths for the singlet and the triplet. Additionally, Aj can be used to tune the coupling 
strength of the triplet independently of the singlet. According to equations (S-57)-(S-23), the 
couplings Aa and Aj also induce vertical polaronic shifts (i.e., of the energy at the minimum) 
which are given by 

A|_ = 4(Ae - Aa)2 , (S-26) 

= (2Ae - IXj f . (S-27) 

Here it is important to note that the triplet polaronic shift does not depend on Aa, and in 
consequence its coupling to the vibration can be modified through Aa in equation (S-25), whereas 
the polaronic shift in equation (S-27) remains unaffected. The coupling Aj can therefore be used 
to tune the polaronic shift to the triplet without affecting the singlet S- state. It is thus indeed 
possible with our model to describe the central observation in the experiment. 


3. Spin-orbit interaction effects 

In this section we consider the possible effects on the electronic properties of the many-body 
states of the CNT quantum dot due to a static curvature-enhanced spin-orbit (SO) coupling Aso- 
There are two reasons for this: First, one may wonder how the above scenario is affected by Aso 
general. Second, even when Aso is too small as to produce a significant shift in the line positions 
of the differential conductance, the SO mixing turns out to be crucial to give a nonvanishing 
amplitude to the ground singlet line in Fig. 3c. 
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When considering the nanotube axis oriented along the x-direction as we do here, this SO 
interaction enters through the following Hamiltonian^ 

Hso = AsoCJ’x 0 , (S-28) 

where ax and are Pauli matrices in spin and valley subspaces, respectively. Hence, the spin- 
orbit interaction conserves the valley structure, but it couples different spins, yielding spin-flip 
processes. Written in the BA-basis, 


Hso = Aso ^ , 

T(J 


(S-29) 


the spin-orbit term couples the Kramers doublets \D^) |^+) and \D^) |^+)- Here we used 

the compact notation d = —a and r = —r. 

Naively, one expects that for Aso < ^^TqS- the state-dependent vibrational coupling will 
survive. To again develop some better intuition we discuss as before the adiabatic approximation 
where the mechanical displacement Q is considered as classical parameter. In the full numerical 
calculations, however, we always exactly diagonalize the quantum-dot Hamiltonian model. 

a. One-electron states The adiabiatic Hamiltonian matrix in the N = 1 sector contains a SO 
mixing term between the doublets: 


0 0 Aso \ 

(1) 0 (Q) Aso 0 

0 Aso (Q) 0 ’ 

V Aso 0 0 

where the electronic states are now associated to the following adiabatic potentials 

E^(Q) = E^- fkoAl + + \/2A,)2 . 


(S-30) 


(S-31) 


We therefore expect the SO interaction to couple the doublets at the points of intersection of the 
adiabatic potentials as Q is varied. For Aso < A, the minimum of the ground doublet parabola is 
conserved, and one might expect little change for the line positions in the differential conductance. 
In Fig. S-8a we show the one-electron adiabatic potentials including a small spin-orbit coupling 
Aso = 0.125 A = 0.1 meV and compare them to the bare adiabatic potentials . The important 
conclusion to draw here is that in the parameter regime that we will discuss below in Sec. HC 
there is no significant change around the ground state potential energy minimum. 

b. Two-electron states - state dependent Franck-Condon shifts We now discuss the situation 
in the two-electron charge state. The form of the SO Hamilonian in equation (S-29) shows that 
this operator flips both the orbital and spin projection in one of the two electrons in the quantum 
dot. Therefore, the two states Tq and So with electrons with opposite spins in opposite orbitals 
are not affected by this coupling due to Pauli’s exclusion principle. The other states Sr and 
are indeed mixed, as revealed by the Hamiltonian matrix block: 


/EsAQ) 

Aso 

0 

-Aso 

0 

0 \ 

Aso 

Et^{Q) 

0 

0 

0 

-Aso 

0 

0 

Eto{Q) 

0 

0 

0 

-Aso 

0 

0 

(Q) 

0 

Aso 

0 

0 

0 

0 

EsoiQ) 

0 

0 

-Aso 

0 

Aso 

0 

Es+{Q)J 


(S-32) 


Again, we expect the spin-orbit coupling to be only important around the crossings of the adiabatic 
potentials as Q is varied. The above matrix thus mixes the two-particle states as follows: (i) The 
localized singlets Sr will be repelled by the triplets (ii) These triplet states slightly repel each 
other due to a second order process (since there is no direct matrix element connecting them). 

It now depends on the strength of the spin-orbit coupling relative to the vibrational coupling 
energy (shifts of the potential energy minima xhuo) and the splitting between the “bare” electronic 
states (vertical energies at the minima) whether the SO coupling has a negligible impact on our 
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FIG. S-8. Adiabatic potentials including spin-orbit interaction Aso = 0.125A = 0.1 meV. (a) One-particle 
mixed doublets from equation (S-30). (b) Two-particle potentials calculated from equation (S-32) in the 
subspace mixed by the SO coupling. Adiabatic potentials for Aso = 0 are shown in grey for comparison. 
For large Q we can associate each adiabatic potential to the marked electronic states. In both charge 
sectors, the relevant lower minima are slightly affected by Aso- 


mechanism or not: the SO coupling only has a big effect at crossings of potential energies which 
may be far away from the relevant minima. For the estimated experimental parameters, the 
energy difference between the triplet and the ground singlet is noticeable. Moreover, the adiabatic 
potentials related to these states are strongly shifted, meaning that the SO coupling will not have a 
strong impact on the development of the minima of those if Aso < • This is corroborated 

by Fig. S-8b where the two-electron adiabatic potentials are shown. This is in agreement with 
the measurements. However, even in this case the SO coupling plays a role: it is required for the 
explanation of the peak amplitudes in a magnetic field as we discuss in Sec. IIC 2 b. 

For strong spin-orbit coupling, Aso > the situation is quite different: All involved 

Franck-Condon shifts of the vibrational mode become approximately the same, in clear disagree¬ 
ment with the measurements. Roughly speaking, for strong spin-orbit coupling the states S- and 
To are strongly mixed and the difference in their coupling to the vibration is “averaged out”. 

We conclude two things: (i) The SO interaction (S-28)-(S-29) does not generate a state-dependent 
eleetron-vibration eoupling. (ii) When present and strong, the SO interaction rather tends to 
weaken ik merely renormalizing the vibration frequency. 


B. Master equations - tunneling and relaxation 

In this section we describe the employed method for the calculation of the Coulomb-diamond 
stability diagrams shown in Figures 2f - 2g and Fig. 3c of the main article. 


1. Tunnel processes 

In the stationary limit and for weak couplings to the source and drain leads, the occupations 
probabilities Pa in the dot obey the rate equations 

0 = - WbaPa) , (S-33) 

where Wab = + ^ab represents the probability per unit time for a state transition \b) \a) 

in the quantum dot. Since we will restrict ourselves to lowest order contributions in F, the overall 
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scale of Wab, the above rates coincide with those obtained by Fermi’s Golden’s Rule, namely 

Wab = E WL = E (Ea - E,) , (S-34) 

r rrj 

where ff{x) = [1 + ex.p{r]{x — jXr)/kBT)]~^ is the Fermi distribution function for an electron 
[t] = +1) or a hole (77 = —1) and the sum runs over the reservoirs r = {s, d}. The tunnel current 
Ir that flows out of electrode r is calculated through the standard master equation approach in 
the single-electron tunneling regime (SET) 

Ir = E(^« - Nb)W:,pb, (S-35) 

ab 

where Na = {a\N\a) is the electron number in the quantum-dot state \a). Like the state occupation 
probabilities Pa obtained from equation (S-33), the current thus also depends on the tunneling 
rates F^^. The tunnel rates F^^ are related to the tunnel matrix elements (TMEs) for an 

electron with spin a entering {rj = + 1 ) or leaving (77 = — 1 ) the dot and the density of states in 
the leads p. The latter is assumed to be constant in the model (wide-band limit) and hence 


F 


rr] _ 

ab 




a 


(S-36) 


In equation (S-18) we only accounted for the tunnel matrix elements associated to the pure elec¬ 
tronic states |e), i.e. the many-body eigenstates of equation (S-5) labeled by e = D±, S±, Sq^Tq. 
These now need to be extended to the electron-vibration states |e, u) by adding the Frank-Condon 
overlap of the vibrational wave-functions involved in the tunnel event, i.e. 

rpvar] _ rprar] /q 07^ 

^a,u'^b,u \0 o () 

where now a^b = S±^ and is given by equation (S-18). The Franck-Condon 

coefficient strongly depends on the horizontal shift A between the adiabatic potentials of the 
electronic states. If the dot is initially in the state |6, z/), the probability of a transition to a final 
state |a,z/') shifted in A is modulated by^^^^ 






-aV2/ 





(S-38) 


for F > V (replace z/' z/ for z/ > F) and L'jix) is the associated Laguerre polynomial. Since 
the A-shift strongly attenuates the transition probability between the involved electronic states 
at different regimes of the bias, one of the key aspects of the model is that electronic transitions 
to the triplet state allow a change in the number of vibrational quanta while this is exponentially 
suppressed for transitions to the ground singlet We note that also in the presence of vibrations 
the overall tunneling rate F - entering as an overall factor through equation (S-17), equation (S- 
18) and equation (S-37) - merely sets the scale of the current and is irrelevant to the relative 
magnitude of the different excitations which is of interest here. 

We now proceed with the explicit calculation of the TMEs for the 1 2 transitions. Let us 

consider “charging” transitions, i.e., those that start from the N = 1 charge sector [the “discharg¬ 
ing” transitions involve the TMEs = (T^^^^)*]- Eor a strong intrinsic relaxation (see below), 
the relevant transitions are those which begin from |D_,0), i.e., the ground doublet D- and no 
vibration, z/ = 0. Eor (charging) transitions to a doubly occupied state we use p = 1 and find the 
following amplitudes 


rprCT 


rprCT 


rprcr 


Fup{As.,D_)^tL,^Sr-6 

Fi,p{ATQ,D-)f 

F 


a'a •) 


, (S-39) 

= E„ o{I^To,DVlt+<T'^<Tl^<T’l , (S-40) 

(S-41) 


Here we again used the compact notation a = —a and r = —r. We note that the transition to the 
5'+ singlet from the ground doublet is forbidden in the SET regime since it would involve an orbital 
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flip process which is not present in the tunnel Hamiltonian. The same restriction also applies to 
the S- transitions. However, as discussed in the previous section, a small spin-orbit term 

hybridizes the doublets and localized singlet states, making visible the resonance line associated 
to this last transition (line E in Fig. Id of the main article). 


2. Intrinsic relaxation 

In the measured stability diagrams of the main article, we note that electronic transitions 
starting from an excited state are strongly suppressed. In order to account for this effect in the 
most simple way we model the influence of an environmental bath by allowing for relaxation 
processes. Since the main source of excitations is the electron transport we neglect absorption 
processes due to the bath by setting the bath temperature Tb = 0. Furthermore, we allow in each 
charge sector for the relaxation of any energetically higher-lying state h into any lower-lying state 
a, independently of the spin or the orbital distribution of the involved states. The decay rates are 
assumed to he proportional to the energy difference between these states and is assumed to exceed 
the tunneling relaxation rates: 

- Ea), r,,i{E) = r X {E/ 0.2 meV). (S-42) 

The relaxation rate matrix is thus upper triangular and it is specified in units of T. The 
corresponding rates Wff are added to the golden rule rates of equation (S-34) and the master 
equations are solved using these modified rates. This depends little on the details and has the 
main effect of preventing a very strong nonequilibrium state on the quantum dot, enhancing the 
ground state occupation probability in each charge sector. 


C. Comparison with experiment 

In this section we describe how we proceed in finding a unique parameter regime of the model 
which is able to qualitatively explain the experimental data. The key effect of state-dependent 
vibrational coupling is discussed fully in the zero magnetic field case. The finite magnetic field 
experiment brings in some complications due to spin-dependent tunneling and spin-orbit coupling, 
which are, however, not crucial for the central point of the article. Finally, we discuss the depen¬ 
dence on the parameters and identify the problem that a standard Anderson-Holstein model has 
with explaining the observations. 


1. Zero magnetic field 


a. Vibrations As already outlined in the experimental Sec. IB we calculate the capaci¬ 

tive effect induced by the gates and the source and drain leads from Fig. S-9a. The left and right 
resonance lines (red dashed lines in Fig. S-9a), related to the transition between the N = 1 and 
N = 2 ground states D_ S- , are described by the following relations between the top-gate and 
the source-drain voltages 


Vd = - 


= + 


2a 


tg 


I -|- (q^s O^d) 


Vtg - '^sktg 5 


2a 


tg 


I (Ofg Q^d) 


Vtg = mdVtg 


(S-43) 

(S-44) 


By extracting the slopes of these lines (rUg = — I.I and = 0.9) we obtain the following values 
for the lever arms: cttg = 0.5 and ctg — ctd = — O.I. We next extract the electronic parameters, the 
subband splitting A and the exchange energy J. To this end, we use the formulas [cf. equation (S- 
14), (S-45), and (S-46) below)] for the bias distances jxi = 1.2 meV and 112 = 2.7 meV in energy 
units (|e| = I), measured in Fig. S-9a: 


Ml = Etq — Es_ = 2A — JfA, 
M2 = Esq — Es_ = 2A -]- 3J/4, 


(S-45) 

(S-46) 
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FIG. S-9. Resonance lines fitting for N = 1 ^ 2 transitions, (a) Vibrations switched “off” [zoom-in of 
Fig. S-6b and Fig. 2c in the main article]: The measured fii and voltages are related, respectively, to 
electronic transition energies Etq — Es_ and Esq — Es_ • (b) Vibrations switched “on” [zoom-in of Fig. S- 
6d and Fig. 2d in the main article]: The measured fii here corresponds to Eto — Es_ . In both panels, 
dashed lines are visible transitions described by the model: |T)_,0) |*S'-,0) (red), |D-,0) |Tb,z^) 

with = 0 — 4 (blue), |T)_, 0) ^ \So,iy) with z/ = 0 — 2 (yellow). Yellow lines are on top of blue lines (see 
condition 3 in the text). The jiS-, 0) ^ |T)_, 1) transition is shown in panel (b) in dashed magenta and is 
just discernible in the measured differential conductance. 


We obtain A = {fi 2 + 3/ii)/8 = 0.8 meV and J = /i 2 — /^i = 1-5 meV [cf. equation (S-3)]. 

b. Vibrations In the regime where the CNT quantum dot is coupled to the vibrational 

mode, we aim to account for the following features observed in the experimental data of Fig. S-9b 
(Fig. 2d of the main article): 

1 . The experiment shows that the two-electron ground state is given by the singlet S- and the 
first excited state is the triplet Tq. In addition, their line positions remain almost unchanged 
when turning “on” the vibrational coupling. In our model we thus include the base values 
of J and A, fixed by their fitting to the electronic spectrum in Fig. S-9a. When allow these 
energies to become sensitive to the vibration Q through A a and Aj we must thus impose 
the severe restriction that these couplings do not give rise to a polaronic shift. 

2 . The ground singlet S- transition is very weakly coupled to the vibrational mode, the only 

visible lines of the \D_^-u) |5'_,z/') transitions is the “zero-phonon” one, z/ = z/' = 0. We 

thus require a small value for the horizontal shift between the N = 1 and N = 2 ground state 
potential minima, i.e. | <C 1. In contrast, the triplet transition shows a whole series 

of vibrational sidebands and must therefore require a sizeable coupling, i.e., \Atq,d_ \ — 0.5, 
cf. Fig. 4 of the main article. 

3. Although with vibrational coupling turned “off” the excited singlet Sq transition appears 

in Fig. S-9a, no transitions \D-,0) foi* ^ = 1,2,3,... can be distinguished in the 

stability diagram Fig. S-9b when they are turned “on”. The approximate relation J « 2huo 
suggests these transitions are in fact superimposed with the triplet sidebands: the energy 
difference between Sq and Tq, given by the exchange energy J, happens to be commensurate 
with the vibration energy huj to within the experimental thermal line broadenings. We 
require only J = nhw and determine the best fitting integer n (confirming indeed that 
n = 2). 

In applying these constraints to the model parameters, we used energies obtained from lever arms 
that we extracted independently for Fig. S-9b: in this case we obtain from the slopes (rUg = —1.0, 
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md = 0.85) of the ground state lines, giving the following values: atg = 0.46 and a^ — o^d = —0.08. 

For the vibration energy we use the fiuj = 0.85 meV, the fitted value of the mean level spacing 

AF^vib obtained specifically for the upper panel in Fig. 2h of the main article. 

According to constraint (1) for fixed electronic parameters A and J, we need to adjust the 
energy difference between F^To,o and Es_^o to the measured bias 0.9 mV in Fig. S-9b. This 
implies 

Ml = Etq,o — Es_,o = 2A — — + ha;(A|_ — A|.^). (S-47) 

For the following algebra it is convenient to denote the two free parameters hy x = As_ and 

y = Atq and introduce p = {pi — 2 A J j A) j huj —0.382. The first condition for the couplings 
then reads: 

y = x‘^ — p. (S-48) 

The horizontal shift q = Ato,d_ between the triplet and the doublet can be written diS q = y — x/2 
[cf. equation (S-21)-(S-22)] and hence we have 

y=^+q. (S-49) 

From condition 2 we have \x\ 1 and \q\ ^ 1. Bearing in mind this restriction, we find a unique 

solution for x when requiring both the above two equations to hold simultaneously: 

a: = y ± ^ \Jq‘^ + ^P ■ (S-50) 

Notice here that we cannot take \q\ < ^/^p\/4 0.54 since this would imply no real-valued 

solution. On the other hand, we cannot increase \q\ indefinitely either since otherwise the solution 
for X grows and will violate condition 2 [i.e., the sidebands associated to S- would become visible]. 
For q we pick the value g = 0.6 which gives the best overall agreement with the experiment (number 
of visible triplet sidebands), implying x = —0.04 and y = 0.62. 

Finally, we use the condition 3 to obtain the three vibration couplings. This last requirement 
reads with z = A^^ and equation (S-47) 

nhio ~ Esofi - Eto,o = J - - y ^), (S-51) 

where n is an integer number. By using J = 1.5 meV and fiw = 0.85 meV, we obtain + 
J/fujj 2.15 and this implies that n cannot be larger than 2. As mentioned above, n = 2 yields 
approximately the same line position for So as observed before in the pure electronic regime. 
Solving the above equation with n = 2 for z and using the polaronic shifts definitions in terms of 
the A-parameters (equations (S-21) and (S-22)), we obtain the parameter values used in Figures 2d 
and 3c in the main article. 


= 0.28, Aa = 0.32, Aj = -0.22. (S-52) 

We emphasize that the above procedure essentially determines a unique regime of parameters 
consistent with the experimental results. This concerns their qualitative features: conditions 2 
and 3 only define a range of possible values for the polaronic shifts and therefore small deviations 
of the above obtained values for the A-parameters produce similar results. However, an Anderson- 
Holstein type model (Aa = \j = ^) is eertainly not eonsistent with the measurements, see also 
See. IIC 3. Thus, despite the faet that there are several parameters, the experiment imposes strong 
restrietions, in partieular, limiting the ehoiee of vibrational eouplings. Importantly, the transport 
parameters y, n and ( adjust other aspeets of the transport speetrum but do not generate or affeet 
in an essential way the state-dependent vibrational eoupling. In particular at H = 0 are not that 
important. 


2. Magnetic field spectroscopy 

We now consider the predictions of the above described model for magnetic field evolution of 
the dJ/dV^d-peak intensities shown in Fig. 3 of the main article. 
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a. Zeeman effeet on state-dependent vibrational side bands When applying a static magnetic 
field 5, the many-electron states in the quantum dot experience a Zeeman shift. Since the field 
is perpendicular to the CNT axis the orbital splitting due to B can be neglected^The field 
dependence of the doublet and triplet states with spin projection indices a = {+1,-1} and 
m = {+1,0,—!}, respectively, is 

Edi,v = , (S-53) 

Ets',v = bb “ rnguBB , (S-54) 

where E^J, denote the zero-field eigenenergies (S-23). The singlets energies are obviously inde- 
pendent of B. (Regarding the negligble importance of the singlet 5'+, see the remark in the main 
article’s text and the footnote 28 there.) Since the N = 1 quantum dot ground state (| 1^1,0), 
z/ = 0) evolves with —gii^B/2^ the energy cost for adding one electron to the quantum dot depends 
in all cases on the magnetic field strength, no matter which final state is reached: relative to this 
state, the addition energies for the N = 2 states accessible in the sequential tunneling regime (see 
Sec. IIB) are 


AEs_,u = A4°b + , (S-55) 

= (S-56) 

AE^o^, = A4b + IggsB, (S-57) 

AEs„, = A4b + IggeB. (S-58) 


where AE^l = E^l — E^p_ „ are the zero-field values. The ground-state singlet S- and the excited 
triplet Tq lines thus evolve with opposite slopes as function of B. Since the triplet Tq states have 
the interesting vibrational side bands - both in the experiment and in our calculations - their 
field dependence is plotted such that the triplet peak positions stay fixed, both here in Fig. S-IO 
and in Figures 3a - 3c of the main article. In the theoretical plots this is achieved by replacing 
Vtg ^ htg + \gii^B. 

In Fig. S-IOa we show the magnetic-field evolution of the calculated differential conductance 
along the line Vsd = Ftg + 5 mV for spin-independent tunnel barriers, i.e., ^ = 0 for all B (but 
including SO coupling). What the model accounts for at this level relates to the key observation 
made in the measurements: the triplet Tq state (vertical) show vibrational side bands, whereas the 
ground singlet state S- (left most sloped line) does not. Importantly, the other sloped lines should 
not be confused with vibrational side bands of the ground singlet S-: These are the transitions 
into the m = 0 component of the triplet, Tq , or the exeited singlet Sq (the latter is shifted in about 
2huj with respect to the former, see above in Sec. IIC I) and their vibrational side bands. That 
these lines should be there is clear from the theoretical model and its analysis [cf. equation (S-57) 
and equation (S-58)] but is also borne out by noting the difference between the distance A and 
the following ones B-D [just as in the experimental data in Fig. 3a of the main article]. In the 
experiment these Zeeman split-off triplet states in Fig. S-IOa, as well as the singlet Sq are not 
observed in Fig. 3a of the main article. 

b. Suppression of Zeeman split-off states In order to rationalize the observed absence of Zee- 
man split-off states in the experiment we now include a spin-dependence in the tunnel amplitudes 
through the parameter ( [cf. equation (S-I7)] which we take to be T-dependent: 

1 _ p—gu^B/k^T 

^ 1 -(- f.-gmBjk’SiT ’ (S-59) 

This phenomenological function ensures that for gfi^B ^ k^T we have C = I, i.e., the tunnel 
amplitudes are fully polarized for spin-up carriers, see equation (S-I7), whereas for g/XBB /cbT 
we have C = 0- This spin-dependence is physically motivated by the fact that our quantum dot 
is contacted by CNT leads as shown in Fig. S-3b and explained in the main article [cf. Fig. 2b]. 
When increasing the field by a few Tesla the strong spin-polarization in the tunnel amplitudes 
suppresses the passage of spin-down electrons through the dot and hence the Tq and Sq lines 
(together with their vibrational sidebands) are suppressed since the ground state in the N = 1 
sector is Tt. What remains are the Tq^ excitation and its vibrational sidebands. 



19 



Vtg(mV) Vtg(mV) 


FIG. S-10. Calculated magnetic field evolution for the plotted differential conductance of Fig. 2d of 
the main article: (a) Spin-independent tunnel amplitudes (C = 0 for all B). (b) Spin-dependent tunnel 
amplitudes according to equation (S-59). Further parameters [same as in Fig. 2g of the main article]: 
A = 0.8 meV, J = 1.5 meV, = 0.28, Aa = 0.32, Aj = -0.22, hw = 0.85 meV, Aso = 0.1 meV, = 0.3, 
7 = —0.5 and T = 0.7 K. 


However, the ground singlet S- is also suppressed when assuming this spin-dependence, since 
it requires an additional spin-down electron to fill the bonding orbital (r = —), cf. equation (S-9). 
The presence of the ground singlet S- line, but not its vibrational sidebands is the key feature of 
the experiment. It is at this point that the spin-orbit coupling does have a decisive effect: when 
even a small spin-orbit eoupling is ineluded, the singlet S- reappears by borrowing intensity from 
the triplet Tq (see also Fig. f), but without reinstating the unobserved So and the Zeeman split-off 
states of To and their vibrational sidebands. The SO coupling mostly mixes the ground singlet 
S- and the triplet Tq states [cf. Sec. IIA 3]. In Fig. S-lOb [same data as in Fig. 3c of the main 
article] we show the magnetic field evolution of the calculated differential conductance along the 
line Vsd = Vtg + 5 mV when we include the finite ( of equation (S-59): indeed the suppression 
of the Tq and So lines and their sidebands is maintained while only the ground S- has a clear 
intensity. This produces the observed intensity pattern of Fig. 3a of the main article. It is at this 
last point that the Anderson-Holstein model fails, see Sec. IIC 3. 

From this we infer that in our deviee an interplay of spin-filtering of the suspended CNT parts, 
funetioning as tunnel junetions, and weak spin-orbit interaetion in the CNT quantum dot may be 
responsible for the suppression of Zeeman splitting. As emphasized above it is not responsible for 
the essential effect: The missing vibrational side bands while having a clear singlet ground state 
S- line. 


3. Influence of the various parameters and the problem of Anderson-Holstein coupling 

Having made the detailed comparison with the experiment we now outline the different influence 
the various parameters on the differential conductance. We illustrate this for the case of an 
Anderson-Holstein type model (Aa = Aj = 0) and expose the problem it has in explaining the 
data. 

In Figures S-lla - S-llc we plot the results obtained for our model, keeping all parameters as 
in the main article, except for 7 which we vary and making the Anderson-Holstein approximation 
by setting Aa = Aj = 0 . Also, to obtain a similar number of visible vibrational sidebands we 
have to choose a larger value of A^ = 1 in this case. To find agreement with the experiment the 
lines with positive slope should first of all be suppressed: these are vibrational sidebands relating 
to the \D_,n) ^ I*?-, 0) where n = 1, 2, 3,.. quanta are excited for V = 1. This can be done by 
increasing the junction asymmetry 7 a lot. The best result obtainable this way. Fig. S-llc should 
now be compared with Fig. 2g of the main article: as one adjusts 7 to suppress the vibrational 
sidebands relating to the S- (both with positive and negative slope) the ground state z/ = 0 
transition also becomes suppressed, resulting in a severe current blockade at low bias Vgd ^ 

In contrast, in Fig. 2g these lines are not present due to the state-dependent vibrational coupling 
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FIG. S- 11 . Differential conductance at B = 0 for the Anderson-Holstein model (Aa = Aj = 0 ) with 
Ae = 1. The other parameters are the same as for the state-dependent vibrational coupling model, i.e., 
A = 0.8 meV, J = 1.5 meV, hw = 0.85 meV, k = 0.3, C = 0.0 and T = 0.7 K. The panels show results 
for increasing junction asymmetry 7 : (a) 7 = 0, (b) 7 = —0.5 and (c) 7 = —0.9. For completeness we 
have also included here the SO coupling Aso = 0.1 meV as we also did in our model with state-dependent 
coupling (Aa, Aj 7 ^ 0 ), although this has little effect at H = 0 . 


Aa and Aj which does not block the transport at low voltage. This allows for a much smaller 
asymmetry, leaving the ground-state transition fully visible. 

This problem becomes more pressing when we now consider the magnetic field evolution: while 
S- and its sidebands with negative htg-slope are still present they are quasi-degenerate with the 
triplet excitations at B = 0 in panel (c). Turning on the field B they will however split off from 
the triplet by the Zeeman effect. Also in this case one must assume a spin-dependent tunneling ( 
[equation (S-59)] to find agreement with the experimental data of Fig. 3a where the Zeeman-split 
off states are not observed. As before, this also suppresses the S- transition, both the ground 
one {u = 0 ) and all vibrational sidebands (z/ > 0 ). However, when turning on the SO coupling 
again to remedy this, all S- transition are restored, i.e., including the vibrational side bands 
(this happens again by borrowing intensity from the Tq triplet). This disagrees with the key 
experimental observation that only the triplet transitions show vibrational side bands. 

The above discussion underlines the importance of the experimental advance reported in the 
main article: by being able to switch “on” and “off” the vibrational coupling as well as performing 
magnetic field transport spectroscopy of the quantum dot states, we are able to identify electronic 
states with different coupling to the vibration. 

Next, by varying orbital asymmetry in the tunneling k one changes the magnitude of the S- 
transitions relative to the Tq and Sq transitions: when enhancing k, the tunneling in the antibond¬ 
ing (|—)) orbital is enhanced: since S- has two electrons in that orbital, whereas Tq and Sq have 
only one, this enhances the former relative to the latter. As above, hi cannot be adjusted to find 
agreement with the experiment: as in the discussion of the 7 dependence, it cannot suppress the 
vibrational side bands of S- relative to the S- transition: this requires state-dependent coupling. 

Thus the key problem in trying to use an Anders on-Holstein model to explain state-dependent 
vibration eoupling lies in its basie assumption that all eleetronie states with the same eharge eouple 
equally to the vibration. It fails because it either does not at all show the D_ S- excitation 
(i.e., neither the ground {n = 0 ) or any vibrational sidebands z/ 7 ^ 0 ) or it does show it together 
with all its vibrational sidebands. Which of the two is the case depends on parameters, but an 
apparent state-dependent vibrational coupling seems impossible to achieve. The sidebands for the 
singlet S-, relative to the ground transition, are intense as those of the triplet Tq relative to its 
ground transition. Neither junction ( 7 ), orbital (hi) nor spin-asymmetry {() can get around this 
fact. 
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